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Abstract 

Digeneans are endoparasitic flatworms with complex life cycles including one or two intermediate hosts (first of which is 
always a mollusk) and a vertebrate definitive host. Digeneans may harbor intracellular endosymbiotic bacteria belonging to 
the genus Neorickettsia (order Rickettsiales, family Anaplasmataceae). Some Neoricl<ettsia are able to invade cells of the 
digenean's vertebrate host and are known to cause diseases of wildlife and humans. In this study we report the results of 
screening 771 digenean samples for Neoricl<ettsia collected from various vertebrates in terrestrial, freshwater, brackish, and 
marine habitats in the United States, China and Australia. Neoricl<ettsia were detected using a newly designed real-time PCR 
protocol targeting a 152 bp fragment of the heat shock protein coding gene, GroEL, and verified with nested PCR and 
sequencing of a 1371 bp long region of 16S rRNA. Eight isolates of Neoricl<ettsia have been obtained. Sequence comparison 
and phylogenetic analysis demonstrated that 7 of these isolates, provisionally named Neoricl<ettsia sp. 1-7 (obtained from 
allocreadiid Crepidostomum affine, haploporids Saccocoelioides beauforti and Saccocoelioides lizae, faustulid Bacciger sprenti, 
deropegid Deropegus aspina, a lecithodendriid, and a pleurogenid) represent new genotypes and one (obtained from 
Metagonimoides oregonensis) was identical to a published sequence of Neoricl<ettsia known as SF agent. All digenean species 
reported in this study represent new host records. Three of the 6 digenean families (Haploporidae, Pleurogenidae, and 
Faustulidae) are also reported for the first time as hosts of Neoricl<ettsia. We have detected Neoricl<ettsia in digeneans from 
China and Australia for the first time based on PCR and sequencing evidence. Our findings suggest that further surveys from 
broader geographic regions and wider selection of digenean taxa are likely to reveal new Neoricl<ettsia lineages as well as 
new digenean host associations. 
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Introduction 

Neorickettsia (family Anaplasmataceae) is a genus containing 
obligate intracellular endosymbionts of digeneans (Platyhelmin- 
thes, Digenea). Although relatively small, this genus has received 
significant attention in the recent years, with a rapid increase in 
the number of known species level lineages in this group [1], [2], 
[3], [4], [5]. Meorickettsia are vertically transmitted through all 
stages of complex digenean life cycles. Additionally, they are 
capable of being horizontally transmitted to the vertebrate hosts of 
the digenean, both human and wildlife, where they can cause 
disease [1], [4], [6], [7], [8], [9]. These diseases are potentially 
debilitating, e.g., Sennetsu fever in humans (Neorickettsia sennetsu), or 
even fatal, e.g., salmon dog poisoning [Neorickettsia helminthoeca), and 
Potomac horse fever [Neorickettsia risticiij [4]. 

Neorickettsia have now been reported from several countries and 
all continents including Antarctica [4], [5] . However, a majority of 
the records are based almost exclusively on immunological 



detection of Neorickettsia in vertebrate hosts, mainly horses. It 
should be noted that positive immunological tests of horses could 
be a result of either actual infection or a previous vaccination 
against Neorickettsia risticii (causative agent of Potomac horse fever). 
Therefore, immunological results alone without PCR-based 
confirmation need to be considered with caution [10], [11]. At 
the same time, screening of digeneans for Neorickettsia has been 
limited, which prevents understanding of the actual diversity of 
these bacteria and potential sources of infection of vertebrate 
animals including humans. Only 21-22 digenean species, in most 
cases identified only to genus or family level, have been previously 
confirmed as hosts of Neorickettsia [3], [4], [12], [13], [14], [15], 
Records resulting from PCR detection of neorickettsia! DNA are 
known only from North and South America, eastern Asia, and 
Antarctica. However, most of these records are based on 
Neorickettsia detected in vertebrate host tissues (horses, dogs, 
humans, fish). As posited by Vaughan et al. [4], future studies 
focusing on screening digenean extracts for Neorickettsia will reveal 
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additional host associations of these bacteria and new pathways of 
their circulation in nature. This has already been demonstrated in 
a recent study by Tkach et al. [3] who found 4 species level genetic 
lineages of .Neorickettsia in 7 species of digeneans belonging to 7 
different families. 

Currendy, there are 3 named species and 10 not formally 
named genotypes of Neorickettsia that are likely to represent 
additional species based on levels of 16S rRNA sequence 
divergence [2], [3], [4], [5], [16], [17], [18], [19]. Of these 13 
known species/genotypes of Neorickettsia, 6 have been found in 
North America. Circulation of digenean hosts of known Neorick- 
ettsia occurs in either freshwater (Rainbow trout agent, Catfish 
agent 1, Catfish agent 2), freshwater/ terrestrial (JV. helminthoeca, N. 
ruticii, Elokomin Fluke Fever agent, Diplostomum agent), or fuUy 
terrestrial [N. risticii) environments [3], [4], [12], [13], [20], [21], 
[22]. Presently, however, there are no records Neorickettsia from 
digeneans having fuUy marine life cycles. Presence ot Neorickettsia in 
tissues of marine fishes in Antarctica and the Gulf of Mexico [5], 
[23] may suggest a marine circulation pathway, however, these 
authors did not screen any digeneans. 

In this study, we used real-time PCR-based detection methods 
to study the diversity of Neorickettsia in previously understudied or 
unstudied regions. To accomplish this goal, we screened for 
Neorickettsia numerous DNA extracts from an extensive collection of 
adult digeneans collected from various vertebrates as well as snail 
hosts in the United States (several states), Australia, and China. As 
a result, we have found 7 new genetic lineages oi Neorickettsia some 
of which may potentially represent new species. AH of our records 
represent new digenean host associations. Our findings also 
expand the range of circulation pathways known for Neorickettsia. 
We have for the first time detected Neorickettsia from Australia and 
China using PCR-based detection and DNA sequencing. We have 
also conducted a molecular phylogenetic analysis in order to 
estimate interrelationships among the newly discovered genotypes 
with previously known named species and not yet named lineages 
of jV«oncfetoa. 

Materials and Methods 

Ethics statement 

Digeneans were collected from marine and freshwater fishes, 
amphibians, reptiles, birds, mammals, and invertebrates from 
2009 through 2012 from multiple localities in China, Australia, 
Argentina, Costa Rica, and several states in the USA (Mississippi, 
Louisiana, Florida, Oregon, North Deikota, Minnesota). An 
lACUC protocol 10100105 was issued by the University of 
Southern Mississippi for collecting and humane euthanasia of wild 
animals in Mississippi, Florida, and Louisiana. Birds were collected 
by shotgun in accordance to the state and federal permits, 
raccoons were live trapped and euthanized in the field by firearm 
(.22 caliber rifle) in accordance to issued permits, and amphibians 
and reptiles were collected by hand and humanly euthanized by 
immersion in a solution of chlorobutanol (chloretone) in water in 
accordance with the issued permits and lACUC protocol. lACUC 
protocols were not required for invertebrate collecting in Oregon, 
North Dakota, and Minnesota or for fish. Fish were purchased 
dead from food markets in Argentina, Costa Rica, and China. In 
Australia, fish were caught by cast net by one of the authors (EP) 
from around Townsville (Queensland), Darwin (Northern Terri- 
tory) and Broome (Western Australia). Fish w("re placed on ice- 
after capture, and if necessary were sacrificed with the use of a 
250 mg/1 bath of Tricaine methane sulfonate (MS-222) in 
accordance with the issued permits. 



Sample collections 

AU digenean samples were collected over a period of 7 years as 
parts of several independent projects dealing with parasite 
biodiversity and systematics. Thus, the source of DNA available 
for Neorickettsia screening was opportunistic, but as inclusive as 
possible. Required scientific collecting permits were obtained in all 
cases. The Florida Fish and Wildlife Conservation Commission 
issued die permits LSSC-1 1-00074 and FNW-13-05(renewal) that 
allowed for the collection of invertebrates, amphibians, reptiles, 
and mammals as well as freshwater fishes in Florida. The 
Mississippi Department of Wildlife Fisheries and Parks issued a 
permit 0123131 for the collection of invertebrates, amphibians, 
reptiles, mammals, and freshwater and marine fishes in Mis- 
sissippi, and the Louisiana Wildlife and Fisheries issued a permit 
LNHP-13-017 also for the collection of the above mentioned 
animals in Louisiana. Additionally, a federal permit (MB68 1207-0) 
was issued by the US Fish and wildlife service for collection of 
reptiles, mammals, birds and freshwater fish in Mississippi and 
Louisiana. Three separate permits were issued in Australia for the 
collection of fish; Western Australia (Government of Western 
Australia, Department of Fisheries, SPA-01-10); Northern Terri- 
tory (Northern Territory' Government, Department of Resources, 
SI 7/2932); and Queensland (Department of Primary Industries 
and Fisheries 133621). Only abundant commercial fish species 
were purchased from local markets and fishermen in China, 
Argentina, and Costa Rica, and therefore, permits were not 
required. Additionally, aquatic snails were collected from North 
Dakota, Minnesota, and Oregon. Snails were identified using 
Burch's "North American Freshwater snails" [24], and also did not 
require permits. We do not provide a complete list of screened 
digenean species because in the case Neorickettsia a negative result 
does not necessarily mean that a certain species of digenean 
cannot be a host for Neorickettsia. As discussed by Tkach et al. [3] 
digeneans do not have any known co-dependency with their 
Neorickettsia endosymbionts and therefore, different individuals of 
the same digenean species may or may not be infected. 

Sample processing 

Live worms were rinsed in saline, examined briefly, killed with 
hot water, and fixed in 70'!^) ethanol that allowed for both 
morphological examination and molecular study. When numerous 
specimens were found, some were fixed in 95% ethanol for 
molecular work. Fixed worms were stained in aqueous alum 
carmine, Mayer's hematoxylin, or Van Cleave's hematoxylin; 
dehydrated in a graded ethanol series; cleared in clove oil (carmine 
and Van Cleave's) or methyl salicylate (Mayer's); and mounted 
permanently in Damar Gum for morphological identification. 

Genomic DNA was extracted from individual adult worms and 
pooled larval stages either using the Qiagen DNAeasy tissue kit 
(Qiagen, Inc., Valencia, California) following the manufacturer's 
instructions or the guanidine thiocyanate method according to 
Tkach and Pawlowski [25]. 

Molecular screening 

DNA extracts were first tested for the presence of Neorickettsia 
using a real-time PGR protocol designed and extensively tested by 
one of the authors (SEG). Five microliters of each DNA extract 
were used. The real-time PGR amplified a 152-bp portion of the 
3 ' end of the heat shock protein coding gene, GroEL. The primer 
pair (designed by SEG) used is listed in Table 1. Samples that 
tested positive with real-time PGR were verified using a 
substantially modified nested PGR protocol initially described by 
Barlough et al.[26]. Five microliters of each DNA extract were 
used for the first PGR reaction and 1 (tl of the first PGR product 
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was used for the nested PGR. A 1470 bp long fragment of the 16S 
rRNA gene was first amplified using the primer pair listed in 
Table 1, designed by SEG. The nested PGR step amplified a 1371- 
bp fragment using internal primers also listed in Table 1. The 
same nested PGR primers were used in sequencing reactions along 
with internal forward and internal reverse primers (Table 1), 
designed by SEG. 

The real-time PGR reactions were run on a Bio-Rad GFX96 
Touch real-time PGR detection system (Bio-Rad Laboratories, 
Hercules, GA) using iTaq universal sybr green supermix (Bio-Rad 
Laboratories, Hercules, GA) according to the manufacturer's 
instructions. A two-step program was used with a denaturation 
temperature of 95°G for 3 seconds, followed by an annealing and 
extension temperature of 64°G for 25 seconds and 36 cycles. In 
addition, a melt curve was run starting from 72°G and increasing 
at 0.2°C increments every 5 seconds until reaching 87°G. The 
nested PGR reactions were run on EP Gradient thermocycler 
(Eppendorf, Hauppauge, NY) using Quick load OneTaq master- 
mbc (New England Biolabs, Ipswich, MA) according to the 
manufacturer's instructions. Annealing temperature of 56°G and 
40 cycles were used in both first and nested PGRs. DNA of jV. 
semetsu used as a positive control was graciously provided by Dr. 
Sabine Dittrich (Lao Oxford Mahosot Wellcome Trust Research 
Unit). Pure water was used for negative controls in both real-time 
and nested PGRs. 

Digenean host identification 

Adult digeneans were identified based on their morphology 
using stained whole mounts. Identification to generic level was 
done using the three-volume "Keys to the Trematoda" 
[27] [28] [29] with subsequent species identification done using 
numerous original descriptions and revision articles. Meorickettsia 
positive cercariae or metacercariae were identified to the lowest 
possible taxonomic level using a partial sequence of the nuclear 
large ribosomal subunit gene (28S). Digenean DNA was amplified 
from the extracted DNA by PGR using the primer pair listed in 
Table 1. The same PGR primers and additional internal primers 
were used for sequencing (Table 1). 



DNA sequencing 

PGR amplicons of both J^eorickettsia and digeneans were purified 
using the Zymo DNA Glean & Goncentrator -5(Zymo Research, 
Irvine, GA) or ExoSap PGR clean-up enzymatic kit from 
Affimetrix (Santa Glara, GA) according to the manufacturer's 
instructions. The PGR products were cycle-sequenced using ABI 
BigDye chemistry, ethanol precipitated, and run on an ABI Prism 
3100 automated capillary sequencer. Gontiguous sequences of 
Neorickettsia and digeneans were assembled using Sequencher ver. 
4.2 (GeneGodes Gorp., Ann Arbor, MI) and submitted to 
GenBank under accession numbers KF661342-KF661351 and 
KF878083- KF878084. 

Phylogenetic analysis 

Newly obtained sequences (Table 2 and Fig. 1) and sequences of 
neorickettsiae from GenBank were used in the phylogenetic 
analysis. The new sequences and sequences obtained from 
GenBank were initially aligned with the aid of GlustalW as 
implemented in the BioEdit program, version 7.0.1 [30]. The 
alignments were manually refined in MacGlade, version 4 [31]. 

Phylogenetic analysis was carried out using Bayesian inference 
(BI) as implemented in the MrBayes program, version 2.01 [32] 
with the following nucleotide substitution parameters: Iset nst = 6, 
rates = invgamma, ngammacat = 4, that correspond to a general 
time reversible (GTR) model including estimates of the proportion 
of invariant sites (I) and gamma (G) distributed among-site rate 
variation. Posterior probabilities were approximated over 
1,500,000 generations for both data-sets, log-likelUiood scores 
plotted and only the final 75% of trees were used to produce the 
consensus trees by setting the "burnin" parameters at 375,000 
generations. The GTR+I+G model was used for the analysis based 
olf of the results obtained from JModelTest, version 0.1.1 [33], 
[34]. 

Results 

A total of 7 75 digenean samples were screened for Neorickettsia. A 
list of digenean families screened and the number of extracts 
corresponding to each family are provided in Table 3. Screening 
revealed 8 different genetic lineages ot Neorickettsia in 8 digenean 



Table 1. PGR Primers used in the study. 



Reaction type 


Primer 


Sequence (5'-3') 


Real-time 


groel-lSOOF 


ATAGATCCAGCKAAGGTAGTGCGTGT 




groel-1620R 


TTCCACCCATGCCACCACCAGGCATCATTG 


l" round PGR {Neorickettsia) 


nl6S-25F 


TCAGAACGAACGCTAGCGGT 




nl6S-1500R 


AAAGGAGGTAATCCAGCCGCAGGTTCAC 


Nested PCR/sequencing {Neoricl<ettsia) 


nl6s-50F 


TAGGCTTAACACATGCAAGTCGAACG 




nl6S-1400R 


CGGTTAGCTCACTAGCTTCGAGTAA 


Sequencing {Neoricl<ettsia) 


16S-n900F 


GACTCGCACAAGCGGTGGAGTAT 




16S-n900R 


ATACTCCACCGCTTGTGCGAGTC 


PCR/sequencing (digenean) 


digl2 


AAGCATATCACTAAGCGG 




1500R 


GCTATCCTGAGGGAAACTTCG 


Sequencing (digenean) 


300F 


CAAGTACCGTGAGGGAAAGTTG 




900F 


CCGTCTTGAAACACGGACCAAG 




300R 


CAACTTTCCCTCACGGTACTTG 




ECD2 


CTTGGTCCGTGTTTCAAGACGGG 



doi:1 0.1 371 /journa!.pone.0098453.t001 
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■ Wolbachia pipientis (JX1 82384) 



•Anaplasmataceaesp.(Cstriata)(EU780453) 

'Neorickettsia sp. (Rainbow trout agent) (AF206298) 
•Neorickettsia sp. (Catfish agent 2) (KF661 349) 

■ Neorickettsia sp. 1 1Crepidostomum affine) (KF661 342)* 

■ Neorickettsia sp. 2 [Saccocoelioides beaufbrti) (KF661 343)* 
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Neorickettsia sp. (/M. armatus) (EU780452) 
Neorickettsia sp. (W. coriiceps) (FJ456799) 
00 ^Neorickettsia sp. {N. coriiceps) (FJ456797) 
I — Neorickettsia sp. ICaceratus) (FJ456595) 
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lOOfNeorickettsia risticii (Illinois strain) (NC01 3009) 
^Neorickettsia risticii (AF380257) 

[Neoricitettsia sp. (Metagonimoides sp.) (KF661 347)* 
Neorickettsia sp. (SF agent) (U34280) 

{Neorickettsia sp. 3 (Pleurogenidae) (KF661 344)* 
/Veor;ckettsfasp.6 (Lecithodendriidae) (KF878083)* 
-A/eorfctetts/ar/st;c/i'(Caddisfly larvae) (AF206300) 
>- Neorickettsia sp. 7 {Deropegus aspina) (KF878084)* 
— Neorickettsia sennetsu (M73225) 
Neorickettsia risticii (Plagiorctiis elegans) (KF661 351 ) 
— Neorickettsia sp. 4 (Saccocoelioides lizae) (KF661 345)* 
• Neorickettsia sp. 5 [Bacciger sprenti) (KF661 346)* 



87 



B 

90 



Figure 1. Phylogenetic relationships among 25 taxa of bacterial endosymbionts in the family Anaplasmataceae resulting from 
Bayesian analysis. Phylogenetic relationships among 25 taxa of bacterial endosymbionts in the family Anaplasmataceae resulting from Bayesian 
analysis (1,500,000 generations) of partial sequences of 16S rDNA gene. Posterior probabilities greater than 80% are shown above internodes. Roman 
numerals (I, II, II) represent the different clades within the "Neorickettsia clade" and letters (A, B, C) correspond to the subclades within clade III. An 
asterisk (*) at the end of a taxon name indicates a new genotype of Neorickettsia discovered in this study. GenBank numbers are given here for all 
taxa. 

doi:1 0.1 371/journal.pone.0098453.g001 



species belonging to 7 different families (Table 2). Neorickettsia 
infections were detected in the foUovving collecting sites: Pearl 
River, Mississippi (30°20'36"N, 89°38'O3"\\0; Cedar Key, Florida 
(29°08'13"N, 83°02'36"W); Williford Springs, Florida 
(30°26'20.40"N, 85°32'50.77"W and 30° 26'21.95"N, 
85°32'42.26"W); Daya Bay, Guangdong, China (22°43'N 
114°32'E); EH Creek, Queensland, Australia (25°15'45"S, 
152°48'28"E). We have detected Neorickettsia from digeneans for 
the first time from the Australian continent and China as well as 
from Orc-gon and Florida. One of the forms discovered in our 
study is identical to the SF agent (Table 2). Sequences of the 7 
other forms clearly differed from all other known forms of 
Neorickettsia, therefore, potentially representing new species. 

The phylogentic analysis was run using a dataset including 25 
ingroup taxa and Wolhachia pipientis (a member of Anaplasmata- 
ceae endosymbiotic in insects and fUariid nematodes) as an 
outgroup. The alignment included a total of 1,258 sites, of which 
1,249 could be aligned unambiguously. Positions that could not be 
aligned unambiguously were excluded from the analysis. Bayesian 
analysis produced a tree where all Neorickettsia sequences clearly fall 
into a well-defined clade, with a 100% support (Figs. 1 and 2). 



The "Neorickettsia clade" includes three major, well supported 
clades (Figs. 1 and 2). Clade I is represented by a single sequence 
from Bacciger sprenti in Australia. Clade II includes 6 sequences, all 
of them representing yet unnamed lineages/species of Neorickettsia 
found in either digeneans parasitic in fish or in fish tissues. Two of 
these forms are represented by our newly sequenced genotypes, 
Neorickettsia sp. 1 (from Crepidostomum affine in Mississippi) and 
Neorickettsia sp. 2 (from Saccocoeliodes heaiiforti from Florida). All 
internal sub-clades of clade II are well resolved and lOO'/o 
supported. Five lineages of this group are from North America 
while one is from Southeast Asia. 

The major clade III contains a majority of the Neorickettsia 
species, including all three currendy recognized named species of 
these bacteria, namely N. helminthoeca, N. risticii, and N.sennetsu 
(Figs. 1 and 2). This large clade is split into three sub-clades 
indicated as A, B, C on Figures 1, 2. Sub-clade A contains a 
diverse assemblage of lineages that includes N. helminthoeca, 
Neorickettsia sp. from a Diplostomum sp. published by Tkach et al. 
[3] as separate branches, and a well resolved cluster containing 
one form from Southeast Asia and three genotypes from the 
Antarctic [5] . All digenean hosts of the members of this sub-clade 
are parasitic in fish at least at some phase of their life cycle. Sub- 
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Table 3. List of digenean families screened for Neorickettsia and number of individual extracts associated with each family. 



Di^enedn Fsniily 


Numbsi' of GxtfBcts 


Digenedn Fsinily 


Number of extrscts 


Dic|6n63n Fsmily 


NurnbGr of extrscts 


Acanthocolpidae 


2 


Fellodistomoidae 


3 


Paramphistomidae 


2 


Allocreadiidae 


40 


Gorgoderidae 


33 


Philophthalmidae 


5 


Apocreadiidae 


46 


Haploporidae 


227 


Pleurogenidae 


4 


Atracotrematidae 


4 


Haplosplanchnidae 


17 


Pronocephalidae 


4 


Azygiidae 


13 


Hemiuridae 


6 


Psilostomidae 


4 


Bivesculidae 


8 


Heterophyidae 


5 


Renicolidae 


2 


Bucephalidae 


54 


Lecithodendriidae 


17 


Rhopaliidae 


2 


Clinostomatidae 


3 


Lepocreadiidae 


43 


Spirorchiidae 


8 


Cryptogonimidae 


20 


Lissorchiidae 


8 


Telorchiidae 


4 


Cyathocotylidae 


2 


Macroderoididae 


23 


Troglotrematidae 


8 


Derogenidae 


6 


Microphallidae 


28 


Zoogonidae 


3 


Diplostomatidae 


16 


Monorchiidae 


17 


Unidentified 


30 


Echinostomatidae 


1 


Notocotylidae 


1 






Eucotylidae 


1 


Opecoelidae 


22 






Faustulidae 


22 


Opisthorchiidae 


7 







doi:1 0.1 371 /journal.pone.0098453.t003 




Wolbachia pipientis 



- Anaplasmataceae sp. (C. striata) 


OR 


? 


- Neorickettsia sp. (Trout agent) 


NE 




- Neorickettsia sp. (Catfish agent 2) 


NE 




- Neorickettsia sp. 1 {Crepidostomum affine)* 


NE 




- Neorickettsia sp. 2 (Saccocoelioides beauforti)* 


NE 




- Neorickettsia sp. (Catfish agent 1 ) 


NE 




- Neorickettsia sp. (X. cancila) 


OR 


? 


- Neorickettsia helminthoeca 


NE 




- Neorickettsia sp. {Diplostomum sp.) 


NE 




- Neorickettsia sp. {M. armatus) 


OR 


? 


- Neorickettsia sp. (/V. coriiceps) 


AN 


? 


- Neorickettsia sp. {N. coriiceps) 


AN 


? 


- Neorickettsia sp. (C.aceratus) 


AN 


? 


- Neorickettsia risticii (Illinois strain) 


NE 


? 
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Figure 2. Phylogenetic tree resulting from Bayesian analysis (1,500,000 generations) of partial sequences of 16S rDNA gene. 

Phylogenetic tree resulting from Bayesian analysis (1 ,500,000 generations) of partial sequences of 1 6S rDNA gene. Posterior probabilities greater than 
80% are shown. Roman numerals (I, II, II) represent the different clades within the "Neorickettsia clade" and letters (A, B, C) correspond to the 
subclades within clade 111. An asterisk (*) at the end of a taxon name indicates a new genotype of Neorickettsia discovered in this study. Digenean 
vertebrate definitive host groups (fishes, mammals and birds) are indicated by symbols. Zoogeographic regions are shown as the following; NE 
(nearctic), EP (eastern palearctic), OR (oriental), AU (australian), and AN (antarctic). 
doi:1 0.1 371/journal.pone.0098453.g002 
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clade B is 100% supported, however, internal interrelationships of 
its constituent taxa are less resolved than most of the topologies 
elsewhere in the tree. This sub-clade incorporates 4 Neorickettsia 
forms discovered in this study. The basal taxon in this sub-clade is 
the form reported as Neorickettsia risticii by Tkach et al. [3]. The 
closest derived taxon to it is jV. sennetsu, the agent causing human 
disease in Southeast Asia. The remaining lineages of the sub-clade 
B form a polytomy that includes several lineages of jV". risticii and 
the SF agent. Sub-clade C contains only our newly sequenced 
form from China. AU species/ genotypes ot Neorickettsia in this clade 
are associated with digeneans who use mammals and/ or birds as 
their definitive hosts. 

Discussion 

The newly developed real-time PGR protocol has provided for a 
quick and sensitive way of screening large collections of digenean 
DNA (extracts for Neorickettsia. Removing the necessity of running a 
nested PGR and screening the results with gel electrophoreses for 
every extract has greatly reduced the amount of time required to 
screen large collections of DNA extracts. The elimination of the 
use of post-PGR handling steps in real-time PGR also reduces the 
likelihood of false positive results caused by contamination. 
Although real-time PGR using TaqMan probes has been used in 
the detection of jV. risticii trom snails and horses [35], the technique 
was not used for a broad screening of digeneans for Neorickettsia and 
there was no published protocol of "regular" real-time PGR using 
sybr green or another dye as an alternative to the TaqMan probes. 
Another advantage of the real-time PGR protocol is the ability to 
run melting curve analysis, which is particularly useful for 
differentiating potential non-specific binding of primers. It was 
of great importance for our study since we used relatively generic 
primers targeting all Neorickettsia species and genotypes. The 
combination of real-time PGR, nested PGR, and sequencing 
allows for a wry robust and accurate screening. 

Our screening has revealed 8 different genotypes oi Neorickettsia: 
a new genotype 1 from Mississippi, new genotypes 2, 3, and SF 
agent from Florida, new genotypes 6 and 7 from Oregon, a new 
genotype 4 from Ghina, and a new genotype 5 from Australia 
(Table 2). Our findings represent the first records oi Neorickettsia in 
digeneans in Australia, Ghina, and two states (Oregon and Florida) 
in the USA. All 8 digenean species detected as hosts oi Neorickettsia 
in our study (Table 2) represent new host associations for 
Neorickettsia, although Pusterla et al. [13] detected Neorickettsia 
within a species oi Deropegus from Galifornia, neither the species of 
Deropegus or Neoriciiettsia were identified or sequenced. Members of 
the families Haploporidae, Pleurogenidae, and Faustulidae have 
not been previously reported as hosts oi Neorickettsia. 

Sampling of digenean DNA extracts for Neorickettsia was mosdy 
opportunistic. As stated in the Materials and Methods, digeneans 
were collected for various projects prior to our study and were not 
initially intended for detecting Neoriekettsia. This explains the 
unevenness in the distribution of screened samples among 
digenean families. Digeneans were not targeted based on their 
likelihood of harboring the bacterial endosymbionts because very 
litfle is currendy known about the evolutionary associations among 
Neoriekettsia and their digenean hosts. Therefore, we screened as 
many samples and as diverse digenean taxa as possible. Our study 
is the first geographically and taxonomically broad screening of 
digeneans for the presence Neorickettsia. Its purpose is to provide 
baseline data for future more focused studies within certain 
digenean families or those parasitic in certain hosts. 

Our study has revealed a new circulation pathway oi Neorickettsia 
in the natural environment. Until now, Neorickettsia have been 



found in digeneans with entirely freshwater, freshwater/ terrestrial 
or entirely terrestrial fife cycles [3], [4], but not in digeneans with 
completely marine life cycles. Reports ot Neorickettsia in tissues of 
notothenioid fishes (Notothenia coriiceps and Chaenocep/ialus aceratus) in 
Antarctica [5] and mullet Mugil cephalus in the Gulf of Mexico [23] 
may suggest a marine circulation pathway, however, neither of 
these studies screened any digeneans. We have found three 
difierent forms of Neorickettsia (genotypes 2, 4, and 5) from three 
digenean species with completely marine life cycles, namely 
Saccocoeliodes beauforti, Saccocoelioides lizae, and Bacciger sprenti. This 
suggests that Neorickettsia can be encountered in almost every type 
of environment suitable for digenean life cycles. 

One of our Neorickettsia genotypes obtained from Metagonimoides 
oregonensis in Florida fuUy matched the sequence of SF agent. SF 
agent was initially found in metacercariae of the heterophyid 
Stellantchasmus falcatus [36] infecting grey mullet in Japan [13], [37]. 
Our finding is the first record of SF agent outside Japan. Both S. 
falcatus and M. oregonensis belong to the family Heterophyidae, 
therefore, there is a probability that SF agent could be found in 
other members of this large digenean family, members of which 
are parasitic primarily in fish-eating birds and mammals including 
humans. SF agent causes mild clinical symptoms in experimentally 
infected mice and dogs, but not monkeys or humans [8], [11], 
[38], [39], [40], [41]. Thus, it cannot be excluded that it may 
cause yet unknown disease in wildlife. 

Our phylogenetic analysis has produced a tree with strong 
branch support for most topologies (Figs. 1 and 2). Our data 
corroborate the conclusion by Seng ct al. [2] who suggested that 
their sequence obtained from Channa striata most likely represents a 
separate genus of Anaplasmatac(-ac. On the other hand, the 
Anaplasmataceae sp. from Mastacembelus armatus considered by 
Seng et al. [2] to be another new genus, clearly falls into the clade 
II of the "Neoriekettsia clade"(Figs 1, 2). NeorickMsia sp. ,5 genotype is 
of particular interest because it represents a separate clade of 
Neorickettsia with unresolved affinities to other members of the 
genus. Denser sampling of digeneans from both marine and 
freshwater habitats in Australia is necessary to see whether 
Neorickettsia sp. 5 is a unique, divergent form or a member of a 
larger lineage that includes additional, yet undiscovered, taxa. An 
additional systematically rele\'ant result of the present phyl()g(;netic 
analysis is the position of the N. risticii strain from Plagiorchis elegans. 
This form was initially identified by Tkach et al. [3] as N. risticii 
based on the comparison of much shorter DNA sequences and a 
phylogenetic analysis using those shorter sequences and fewer 
taxa. In the analysis by Tkach et al. [3] this genotype appeared in 
a polytomy that included 4 different isolates of N. risticii and N. 
sennetsu. The present analysis based on longer sequences and 
greater number of taxa places "jV. risticii" from P. ekgans as the 
basal lineage in sub-clade B of clade III. It is separated from the 
remaining taxa in this clade by one of the recognized distinct 
species, N. sennetsu. Thus, this genotype likely represents a different, 
novel species oi Neorickettsia which needs to be better characterized 
and difierentiated from formally named taxa. 

The least resolved part of the tree is the polytomy in sub-clade B 
of clade III that includes several genotypes occupying a derived 
position in relation to N. sennetsu (Figs. 1, 2). The relationships 
between these genotypes are likely to be clarified with the inclusion 
of additional genes in future analyses. The last systematic 
consideration stemming from our analysis is the unresolved 
position of N. helminthoeca in relation to other members of the 
sub-clade A of clade III. Unfortunately, the only 16S sequence of 
N. helminthoeca available in the GenBank is of poor quality and 
contains numerous unresolved/problematic positions. We believe 
that obtaining new high quality sequences of jV. helminthoeca wiU 
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help to clarify its aflinities with the other species/genotypes of 
Meorickettsia. 

Some of the clades in the phylogenetic tree show close 
associations among Meorickettsia and the definitive hosts of 
digeneans in which they were found. For instance, the strongly 
supported clade II notably comprises neorickettsiae obtained from 
digeneans parasitic only in fish, mosdy in North America. This 
clade shows the most clear association with a group of definitive 
hosts of digeneans. Another group that demonstrates a distinct 
pattern of associations with definitive hosts of digeneans is the sub- 
clade B of clade III. Digenean hosts of neorickettsiae in this sub- 
clade are parasitic in either birds or mammals with the exception 
of few cases where the vertebrate hosts of digeneans are unknown 
and Neorickettsia sp. 7 found in Deropegits aspina, a digenean that use 
fish and amphibians as a definitive host (Fig. 2). JV. sennetsu, a 
causative agent of the human disease Sennetsu fever in Southeast 
Asia, is one of the species with yet unknown digenean hosts. 
Nevertheless, the phylogenetic tree topology allows us to 
hypothesize that it should be a digenean that uses either mammals 
or birds as a definitive host. 

Patterns of geographic distribution of most of the digenean 
species harboring Morickettsia in this study are varying and in some 
cases little is known. Crepidostomum qffine was recendy described by 
Tkach et al. [42] from Hiodon tergisus collected from Pearl River 
and Pascagoula River drainages in Mississippi, USA. Based on a 
molecular comparison to other species in the genus by Tkach et al. 
[42] , C. qffine is likely endemic to these river basins in Mississippi. 
Therefore, we hypothesize that the distribution of Neorickettsia sp. 1 
is likely to be limited to that of C. affine. 

Mugil cephalus, the host of two digenean species harboring 
Neorickettsia in this study (Table 2), is a globally distributed marine 
fish. However, these two digenean species, Saccocoelioides beauforti 
and Saccocoelioides lizae, are much more restricted in their 
distribution than their fish host. One of them, Saccocoelioides 
beauforti has only been reported from the southeastern United 
States (North Carolina, Louisiana, Mississippi, and Alabama) [43], 
[44]. The otluT species, Saccocoelioides lizae has only l)een found off 
the cost of South-eastern China [45]. Based on the combined 
distribution data of these digeneans and their fish host, it is likely 
that Neorickettsia sp. 2 may be found in other areas along the coast 
of the southeastern United States and that Neorickettsia sp. 4 may be 
found in other costal txjuntric-s in Southeast Asia. 

Deropegus aspina is only known from the Pacific Northwest of the 
United States in salmonid fishes, including Salmo clarki, S. gairdneri, 
Oncorhynchis kisutch, 0. tshawytscha, and a frog, Rana bqyli, [46]. 
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